Conditions of Low Dimensionality for Strongly Interacting Atoms Under a Transverse 

Trap 
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For a dilute atomic gas in a strong transverse trapping potential, one normally expects that, in 
the ground state, the gas will populate only the lowest transverse level. We show, however, that for 
the strongly interacting gas under a Feshbach resonance, the ground state includes a large fraction 
of atoms in excited levels of the trap, even if the gas is very dilute and the trap is very strong. This 
is because the effective atom-molecule coupling is typically enhanced to many times the trap mode 
spacing by an induced confinement along the untrapped dimension(s). Thus one cannot "freeze out" 
the transverse degrees of freedom except under certain conditions. 



I. INTRODUCTION 

Recently, the system of interacting atomic gas at low 
dimensions has attracted considerable interest as it sup- 
ports a wealth of physics 0, 0, Sllll H H . For ultracold 
atoms, low dimensions are typically achieved by apply- 
ing a transverse optical lattice potential Q . With a one 
(or two) dimensional deep optical lattice potential, the 
atoms trapped at different potential wells basically do 
not interact with each other. For weakly interacting gas, 
one can then freeze out the transverse degrees of free- 
dom by assuming the system to be in the lowest level 
of the trap potential. The atomic gas at each individual 
potential well thus behaves as an effective two (or one) 
dimensional system. 

The condition for achieving effective low- 
dimensionality in an optical lattice becomes more 
involved if one has strongly interacting atomic gas. One 
of the most exciting directions in current atomic physics 
is to study this strongly interacting atomic gas, where 
the effective interaction as measured by the atomic 
scattering length can be tuned in the full range via 
a Feshbach resonance through control of an external 
magnetic field Q. The Feshbach resonance basically 
comes from the coupling of atoms in the open collision 
channel to the molecules in the closed channel. If this 
atom-molecule coupling rate becomes larger than the 
level spacing of the transverse trapping potential, one 
cannot assume a low-dimensional system by freezing 
the transverse mode in the lowest level [lOl [ll|. On 
the other hand, intuitively, one may expect that the 
effective atom-molecule coupling rate always decreases 
when we make the gas more dilute, and the transverse 
level spacing increases when we raise the trapping laser 
intensity. So, for a sufficiently dilute gas under a strong 
transverse trap, one can still get the transverse level 
spacing larger than the atom-molecule coupling rate. 
With this expectation, several recent works have studied 
properties of low-dimensional strongly interacting gas 
by (implicitly ^assuming the transverse mode in the trap 
ground state 0, lal . 

In this work, we want to show that the condition of 
low dimensionality for strongly interacting atoms under 
a transverse trap is more subtle than the above simple 



picture. We will show that even if the gas becomes very 
dilute (in the limiting case, one can just have two atoms 
with the gas density tending to zero), for the ground 
state of the system, we will still have a significant (ac- 
tually, dominant) fraction of atomic population in the 
excited transverse levels under a typical wide Feshbach 
resonance. Furthermore, it is very ineffective to reduce 
the transverse excited fraction by increasing the trapping 
laser intensity. For realistic atoms such as '''"K or ^Li, 
even if the trapping potential is increased to some com- 
pletely impractical level, the transverse excited fraction is 
not yet negligible. The basic reason for this unusual phe- 
nomena is that in the low-dimensional trap, there is an in- 
creased tendency for atoms to pair up spatially also along 
the untrapped dimension(s). This effect lends itself to an 
intuitive understanding of the coupling enhancement in 
the dilute limit. Furthermore, with stronger transverse 
traps, the atoms in the induced pairs will become more 
spatially confined also along the untrapped dimension(s). 
So, although the transverse level spacing increases, the 
effective atom-molecule coupling is also significantly en- 
hanced. As a net effect, the atomic population in the 
transverse excited levels becomes pretty insensitive to the 
magnitude of the trapping potential. Although the result 
here does not preclude the possibility of an effective low- 
dimensional Hamiltonian for strongly interacting atoms 
under the transverse trap, it indeed shows that in gen- 
eral, one can not neglect the atomic population in the 
transverse excited levels, and derivation of an effective 
low-dimensional Hamiltonian would be more subtle and 
tricky than one naively expects [13, [l3 • 

In the next section, we give the general formalism for 
strongly interacting atoms under a transverse trap in the 
dilute gas limit. In that limit, the basic picture is cap- 
tured by two-body physics. For the ground state of the 
system, the atoms form pairs, and interaction between 
the pairs become negligible. There have been a few works 
on description of two-body physics of strongly interacting 
atoms in a trap using the single-channel model Tj. Re- 
cently, there have also been descriptions of the problem 
with a more realistic two-channel model for the Fesh- 
bach resonance (in one-dimensional or three-dimensional 
traps [ill Il3|l. The formalism there, however, neglects 
the atomic background scattering. We extend this for- 



malism to include the background scattering, which be- 
comes necessary when the system is outside of the near- 
resonance region. In Sec. Ill, we present our main cal- 
culation results for the condition of low-dimensionality 
under a transverse one-dimensional or two-dimensional 
traps. We will also give detailed studies of the atomic 
population distributions in the transverse levels and in 
the free dimension(s) for the ^"K and ^Li atoms under 
trapping potentials of various intensities. 



II. GENERAL FORMALISM 

To achieve a d-dimensional {d < 3) atomic gas, we as- 
sume a (3 — d)-dimensional optical lattice applied along 
the transverse direction. The lattice potential barrier is 
high so that the atomic interaction between different po- 
tential wells becomes negligible. With a strong lattice 
potential, the atoms at the bottom of the potential wells 
basically see a transverse harmonic trap. The atomic 
gas in each well can then be modeled as in a (3 — d)- 
dimensional harmonic trap of a characteristic trapping 
frequency uj. The atoms are of mass m and possess 
internal states a — {Tii}- We treat the problem of 
strongly interacting atoms across a Feshbach resonance 
using the standard two-channel model [lj|. The short 
range interaction between closed-channel molecules and 
open-channel atoms is modeled by a delta function. The 
Hamiltonian is expressed as H = Hq + Hi , with 



and 



Hj^gt I d^r (*j*|* 



h.c 



Ub 



d^r^l'l'l^ivl/^ 



(1) 

where \1/ (r) is the atomic field operator, $ (r) is the 
molecular field operator, Vb is the bare detuning, gi, is 
the bare atom-molecule coupling rate, and Ub is the bare 
background atomic scattering rate. The bare parameters 
are related to the physical parameters via the standard 
renormalization relations \i3^ : 



u: 



d^k 1 



r-i = i 



Ur. 



Ub = TUp , gb = Tgp , i/p = z/b + T--^- 



(2) 



parameters gp,Up, Vp are determined from the scatter- 
ing data as Up = Airh^abg/m^ gp = ^ Airh? ^coW\abg\ / m^ 
and i>p — ^co{B — Bq) {hco is the difference in magnetic 
moments between the two channels) |lO| . where we have 
assumed that the s-wave scattering length near resonance 



has the form a, ~ 



abg 



1- 



B~Bo 



with 



a-bg 



as the back- 



ground scattering length, W as the resonance width, and 
Bq as the resonance point. 

We assume the gas to be sufficiently dilute so that we 
need consider only two-body physics within each poten- 
tial well. At a very low temperature, two atoms interact 
and form bound atom-pairs. The interaction between the 
atom pairs is negligible in the limit of a very dilute gas. 
The essential physics is then captured by considering the 
state of two atoms under the above interaction Hamilto- 
nian. For the two-body physics, the center-of-mass de- 
gree of freedom is not influenced by the interaction and 
can be separated under a harmonic potential. So we can 
assume the center-of-mass mode is in the ground state of 
the transverse trap and has zero momentum in the free 
dimension(s) . In this center-of-mass frame, expanding 
the field operators '5 (r) and $ (r) in Eq. (1) in terms 
of harmonic modes in the trapped dimensions and plane 
waves in the untrapped dimensions yields 



Hq — hw 2, emkamkcr^mkcr 



mk 



Vb 



bh 



(3) 



where m indexes trap eigenmodes {mi}, z = 1, . . . , 3 — d, 
and k denotes the wave vector in the untrapped dimen- 
sions {fcj}, j = 1, . . . , d. The operators flmko- and b rep- 
resent the corresponding atomic and molecular modes, 
respectively. As there is only one molecular mode in the 
center-of-mass frame, we drop the index of b. We have 
excluded the center-of-mass energy in the Hamiltonian 
(4). The atom relative energy emk is given (in the unit 
of huj) by 



Cmk 



3-d 



3-d 
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(4) 



where at = \fhjmu) characterizes the trap length scale. 
Likewise, the interaction Hamiltonian Hi in terms of 
these modes has the form 



Hi = 

Ub 



9b 



(3-d)/2 
't ^ mnk 



E ^™n («LikT«n-ki^ + h.c. 



ai'-'^L'^ 



Et t 
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where 



where the subscript p denotes physical parameters, ek — 
h^k'^ /2m, and the integral is taken in three dimensions 
with an explicit energy cutoff Ec imposed on two dimen- 

^^^^/E'c/2^/^7rh^. The physical 



0. Then U- 
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For convenience, we define dimensionless bare parameters where ( (s, x) = limAr^^oo J2 -o (" + ^) 



(N+x)- 
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gbo-t 



■3/2 rr -3 



Ub 



Vb = Vb/hio (7) 



and likewise for dimensionless physical parameters g^, 
Up, and Vp. 

A general two-body state for the atoms and the 
molecule can be expressed as [ll| 



the Hurwitz zeta function and F (x) is the gamma func- 
tion. The above set of equations serve as the basic for- 
malism for determining the state of two atoms in a trans- 
verse trap across a Feshbach resonance. If we take the 
one-dimensional case (d — \) and neglect the background 
scattering (let Up = 0), the above equation for E is re- 
duced to the energy equation in Ref. pU^ . where it is 
derived with a different renormalization method. 



I*) = ( /^^^ + J2 ^™»kaLkT«n-ki I l««c). (8) 

\ mnk / 

The coefficients in this ansatz state (normalized to unity) 
are obtained by solving the Schrodinger equation -ff |^) = 
Efkj\'^), which yields 
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where 
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(12) 
(13) 



Thus, Eq. © determines the eigenenergy E, and Eqs. 
(|10|l and (|ll|l give us the eigenstate as a function of the 
eigenenergy. These equations are expressed in terms of 
the bare parameters, and we need to the use the renor- 
malization relation (2) to transfer them into the ones 
with the physical parameters. One can easily check that 
under the relation (2), Zp{E) = gp {vp — E) = 

-1 / " ^2 X -1 

Zb (E) for any E, and [U^ff (E)] = (Up 



u-E 



re// 



U'^^'' (E) — U^ ^a^hu. The divergence in U^ ^ then 

exactly cancels with the divergence in S {E) . As a net 
result, the above equations ©-dJ retain the same form 
upon renormalization — all the bare parameters are re- 
placed by their physical counterparts, and S {E) and 
95* [E] jdE are replaced by 
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III. TRANSVERSE POPULATION 

DISTRIBUTIONS FOR ^'K AND "^LI AND 

CONDITIONS FOR LOW-DIMENSIONALITY 

In this section, we answer the following basic question: 
For a dilute gas under a strong transverse trap, can we as- 
sume the atoms only populate the lowest transverse level 
at low temperature so that the system becomes effec- 
tively low-dimensional? In the extreme limit of a dilute 
gas, we can consider only two atoms in each potential well 
(the density in the free dimension(s) tends to zero in this 
case). We will show that even in this extreme limit the 
excited fraction is typically still significant (dominant) 
for realistic atoms. If the gas density becomes higher, 
the excited fraction surely cannot decrease. 

To answer this question, we use the above formalism 
for the two-atom state, and calculate the population dis- 
tribution in the transverse levels for '^''K and ^Li atoms, 
which are the relevant species for the current experi- 
ments. We take the scattering parameters VF ~ 8 G, 
Ubg ~ 174 as, ^ico ^ IMfiB for "^^K M and W ~ 300 G, 
abg ~ — 1405 a_B, ^J■co — '^fJ^B for ^Li [l3|. With a typical 
trap frequency w ~ 27r x 62 kHz [J], the physical param- 
eters for the atom-molecule interaction are then given by 
gp = 23 (272), Up = 1.7 (-5.5) for 4°K (^Li). To cal- 
culate the transverse population distribution in the sys- 
tem ground state, we first determine the binding energy 
between the atoms using Eq. Q. The binding energy 
Eb = E — {3 — d) /2 (in the unit of huj), where the latter 
term (3 — d) /2 is the energy for two free atoms (exclud- 
ing the contribution from the center-of-mass mode). As 
the binding energy is of interest by itself, we show Eb 
as a function of the magnetic field detuning B — Bq in 

~~ 40j^ g^^^ 6Lj^ j.j;_ 

For ™K, the binding energy saturates on the 



Fig. 



T^ and [T^bjl with d = 1,2 for 
spectively. 

deep BEG side. This effect comes from the positive back- 
ground scattering length of '"'K, and is related to the 
weakly bound state in the open collision channel. The 
detailed discussion of the background scattering effects 
can be found in Ref. [Tg (see also [13 )■ So, outside of 
the near resonance region, the result here is quite differ- 
ent from the one in Ref. |lj||, where the binding energy 
is calculated for the d = I case without the background 
scattering contribution. From Fig. ^ one can also see 
that with a transverse trap, the binding energy \Eb\ is 
always positive, and only tends to zero when one goes to 
the deep BCS hmit 0,13, 111 (or when B - Bq ^ W ior 
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FIG. 1: Binding energy vs detuning at a; ~ 27r x 62 kHz. Bo is the resonant point in the absence of an optical lattice. The 
inserts show a close-up of the binding energy in the near-resonance region. 
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FIG. 2: Ground state composition vs. detuning. Only the six most significant components are shown explicitly. The label 
(m; n) denotes the sum of Pmn ~ X^k ^mnk ^-nd all components identical by symmetry. 



'^''K, where the scattering length goes to zero). This is 
distinct from the case without traps, where Eb — for 
B>Bo. 

In Fig. 121 we show the population distribution in trans- 



verse levels for '**'K and ^Li in two or one dimensional 
traps as a function of the magnetic field detuning. The 
population fraction in the transverse modes (m; n) is de- 
fined as Pmn = X^k l'7mnk| ■ In Fig. 13 wc show a more 
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FIG. 3: The matrix X^k'^mnk for 2D free atoms on the BEG 
side. 



complete picture of the population distribution at a fixed 
detuning on the BEC side for atoms in a one dimensional 
trap. From the figures, one can see that in general many 
transverse modes are populated. For a fixed mode, the 
population still goes down as the energy of the mode goes 
up, but there are so many excited transverse modes that 
the total population fraction in the excited levels actually 
dominates in typical configurations. 

In Fig. ^ we draw the overall fraction populating 
the transverse excited levels, which is defined as Pex = 
1 — (3"^ — ^j^ ??ook- This fraction needs to satisfy Pex ^ 1 
for the assumption that the atoms only populate the low- 
est transverse level. From the figure, one can see that 
this condition is in general not satisfied for "'"K and ^Li, 
except in the deep BCS limit with the binding energy 
\Eb\ <C 1- For the case of ^Li, if one goes to the very 
deep BEC limit with the closed channel population (the 
bare molecule fraction) (3^1, the condition Pex <C 1 
is also automatically satisfied. For '"'K, because of its 
positive background scattering length, the bare molecule 
fraction /3 remains small even if one goes to the deep BEC 
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FIG. 4; Excited fraction and binding energy vs. detuning for 
cj ~ 27r X 62 kHz. 
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FIG. 5: Excited fraction vs. binding energy for various coupling rates. 



limit [TsL [T^ ■ So the excited fraction Pex continuously 
goes up as one increases the negative detuning. From 
this calculation, it is clear that for the entire region near 
resonance, which is of particular experimental interest, 
one cannot neglect the transverse excited fractions for 
any atoms with a wide Feshbach resonance. 

The condition Pex ^ 1 can only be satisfied in the 
near resonance region for the narrow Feshbach resonance. 
In Fig. |S1 we show the excited fractions for various 
atom-molecule coupling rates. The background scatter- 
ing length still takes the positive or negative values corre- 
sponding to '"^K or ^Li atoms, but we vary the coupling 
rate gp (so the resonance width W changes). One can 
see that the condition P^x <C 1 is satisfied in the whole 
region only when gp < 1. As gp oc VW, the condition 
gp < 1 requires a very narrow resonance with the reso- 
nance width W < O.OIG. On the other hand, one can 
also see from the figure that to satisfy P^x ^ 1, the back- 



ground interaction \Up\ can be somewhat larger than 1, 
but it cannot be arbitrarily larger. We have tested (not 
shown in the figure) that the condition Pex <C 1 breaks 
down when Up is on the order of a few tens. 

The above calculations are done with a fixed trap fre- 
quency w ~ 27r X 62 kHz, as it is typical for current ex- 
periments I3] ■ One may expect that if the trap frequency 
further increases, it will become much easier to suppress 
the transverse excited fraction and to satisfy the condi- 
tion Pex ^1. To test whether this is true, we show in 
Fig. the transverse excited fraction versus the mag- 
netic field detuning for various magnitudes of the trap 
frequency. One can see that even if the trap frequency 
increases by several orders of magnitudes, the transverse 
excited fraction P^x does not change much and remains 
significant. To understand this puzzling effect, we note 
that although the transverse level spacing increases a lot, 
the effective atom-molecule coupling rate also increases. 
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FIG. 6: Excited fraction vs. detuning for various trapping frequencies. 



One cannot assume a constant atom-molecule coupling 
rate as an enhanced trapping in the transverse direc- 
tion will suppress the pair size in all directions including 
the free dimensions (see the next paragraph). The effec- 
tive atom-molecule coupling rate increases significantly 
when the pair size shrinks. As a net effect, the ratio 
between the atom-molecule coupling and the transverse 
level spacing is only a slowly- varying function of the trap 
frequency. The dimensionless parameter g^, which mea- 
sures the effective atom-molecule coupling strength in the 
unit of the trap frequency, actually characterizes this ra- 
tio. From its expression, one can verify that g^ goes as 
w"^/"*. So, the effective ratio g^ only drops by a factor of 
10 (it will still be larger than 1 for ^Li and ""^"K) even if 
the trap frequency increases by four orders of magnitudes 
from its current value (which is almost impossible) . This 
explains the relative insensitivity of the excited fraction 
Pea; to the trap frequency. We can also conclude here 



that for any reasonable trap strength one can imagine, it 
is impossible to neglect the transverse excited fractions 
for realistic atoms such as ^Li and ^°K with a wide Fes- 
hbach resonance. 

With an increased transverse trapping, it is easy to 
understand that the pair size along that direction will be 
suppressed. However, it is not so obvious that the pair 
size along the free dimension(s), where there is no trap, 
will correspondingly shrink. The latter actually comes 
from the interaction effect. Under strong interaction, the 
pair size cannot change only along one direction. To see 
the pair size shrinking along the untrapped directions, we 
can take the Fourier transform of the pair wavefunction 
jyinnk with respect to the momentum k in the free dimen- 
sion(s). This Fourier transform, denoted as yymnr, yields 
the real-space wavefunction, which is the probability am- 
plitude for an atom pair in transverse modes m and n to 
be spaced a distance r apart in free diinension(s). From 
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FIG. 7: Characteristic atomic pair size vs. trapping frequency for a fixed detuning. 



Eq. jni), »ymnr IS given by 
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where Kq {x) — J^ dt cos {xt) / Vt^Tj-l is a modified 
Bessel function of the second kind, and 



^n 



at 



-Eb + J2i=l (™i + "0 



(17) 



which characterizes the pair size in the free diniension(s). 
The dependence of ^mn on the trapping frequency is 
shown in Fig. [T] for different transverse modes (mn) 
on the BEC side. Note that under typical configura- 
tions, the atom population is broadly distributed over 
the transverse modes (as illustrated by Fig. ISJ, so the 
exact pair size in the free dimension(s) should come from 
the average of ^mn over different (mn). From Fig. [7| 
one can see that for typical values of (mn) , the pair size 
significantly shrinks with increase of the trap frequency. 



so its average will follow the same trend. As one moves 
towards the ECS side of resonance, the trend only be- 
comes more dramatic. This provides a physical inter- 
pretation of the coupling enhancement: The transverse 
trap raises the minimum energy of the atoms relative to 
the molecules. Thus, the smaller atom pairs which have 
stronger mixture with the molecules become more ener- 
getically favorable. The net result is a sort of induced 
pair- wise confinement along the free dimension(s), and 
this increased local density causes the effective coupling 
to increase. 



IV. SUMMARY 

In summary, we have shown that for experimentally 
relevant cases, atoms trapped along one or two dimen- 
sions cannot be considered to be in the transverse ground 
level except when well out of the BCS-BEC crossover re- 
gion. In the crossover region and for the ground state 
of the system, a significant fraction of atomic population 



resides in the excited transverse levels, even if the gas is 
very dilute. Furthermore, one cannot effectively suppress 
the transverse excited fraction by raising the trap inten- 
sity. Even with an extremely strong trap far beyond the 
current experimental technology, the transverse excited 
fraction is not negligible yet for realistic atoms such as ^Li 
and '^''K across a wide Feshbach resonance. The conclu- 
sion here is that in the experimentally interesting region, 
one cannot describe strongly interacting atoms under a 
transverse trap as a low-dimensional system by assuming 
a fixed transverse mode. Although this result does not 



exclude an effective low-dimensional description of this 
strongly interacting system, it indeed shows that the ef- 
fective description will become much more subtle, and 
such a description needs to take into account the broad 
population distribution of the atoms in all the transverse 
modes. The derivation of such an effective description 
will be a topic for further investigation |12| . 
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